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Simulation of Impact on A Ductile Polymer Plate 

Rebecca Cremona, Brown University 
Jeffrey A. Hinkley, NASA Langley Research Center 


Abstract 


Explicit finite element calculations were used to visualize the deformation and temperature rise 
in an elastic-plastic plate impacted by a rigid projectile. Results were compared to results of 
experiments involving ballistic penetration of a “self-healing” thermoplastic. The calculated 
temperature rise agreed well with the experimental observation, but the total energy absorbed in 
the penetration event was underestimated in the calculation, which neglected friction. 

Introduction 


When plants and animals sustain damage, they can often heal and regain some of their lost 
function. Inspired by this capability of living systems, technologists have sought to design “self- 
healing” synthetic materials. Several lines of research are being pursued: composites that 
contain crack-filling reagents, self-assembled layers of various sorts, and — the subject of this 
report — thermoplastics that flow to seal a ballistic puncture (1,2). The initial work employed 
Surlyn® 8920 (a duPont™ product). Since then, several other such self-healing thermoplastics 
have been identified (3), but the healing mechanism is not fully understood. Because the healing 
occurs in much less than one second, few experimental techniques are applicable. The present 
report, therefore, applies finite element (FE) analysis in an effort to identify important features of 
the impact event. 

Methods 


TOCHNOG, a free open-source software package (4), was used. Although not as versatile as 
some commercial packages, it has a number of useful features including contact (implemented 
with a penalty formulation), new mesh generation to correct distorted elements, and the ability to 
handle thermal convection and diffusion as well as material velocities and stresses. In 
TOCHNOG, velocities, rather than displacements, are the primary variables and variables 
(assumed to be continuous) are stored in element nodes. An explicit integration scheme was 
chosen. A sample TOCHNOG script is given in the Appendix. 

The polymer constitutive law was elastic-perfectly plastic with a von Mises yield criterion. 
Previous work (5) has shown that this is a good starting point for this type of simulation. Initial 
properties were chosen to be similar to those of Surlyn® at room temperature (Table 1). Models 
were 2-dimensional with axisymmetry. Friction was neglected. Nodes at the circumference of 
the target were fixed to simulate a clamped boundary and the horizontal velocity was constrained 
to zero along the centerline. For most of the calculations, the plate radius was 15mm, i.e. about 3 
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times the radius of the rigid projectile. The initial conditions corresponded to the undeformed 
state (zero strain and temperature). 


T able 1. Initial Assumed Properties of Target material 


Young’s Modulus, E, MPa 

350 

Poisson ’s Ratio, v 

0.4 

Tensile yield stress, cr ys , MPa 

15 

Density, kg/m 2 

950 

Tensile failure strain, % 

470 


Finite element results were visualized with the postprocessor in GiD (6). 
Results and Discussion 


In the experimental ballistic tests described in Reference 1, the only data available are a) the 
change in bullet velocity, b) full-field infrared thermographic temperatures obtained a few tenths 
of a second after penetration, and c) the appearance of the specimen post-impact. Taken 
together, this is not enough information to quantitatively validate a finite element model. 
Nevertheless, it is of interest to study the effects of the model assumptions in a parametric way. 

Dynamics of Impact: As a first demonstration of the software’s capability, consider a model of a 
50 mm-radius x 6 mm-thick plate impacted by a rigid sphere moving downward at a constant 
speed of 2.6 m/sec. After 3.5x10 3 sec, (Figure 1), there is a localized dent, as well as significant 
global bending and stretching of the plate. Due to the axisymmetry, only half the plate is shown. 
When the projectile velocity is 260 m/sec (characteristic of a pistol bullet), on the other hand, 
Figure 2 illustrates the result. At the ballistic rate, the polymer deformation is highly localized 
and the global plate response appears to be negligible. The amount of projectile travel in Figure 
2 is the same as in Figure 1 although the scale of the visualization is slightly different. 

Deformation and Failure: A complete description of the plate penetration would require i) a 
realistic constitutive model and ii) a failure criterion. To simulate healing, the material 
description would have to include in addition Hi) nonlinear viscoelastic recovery, iv) melting, 
v) surface tension, and vi) interdiffusion. While finite element methods are able in principle to 
handle these effects, only the initial deformation is treated in this report. 

Investigators of the self-healing thermoplastics have considered various different projectiles: 

Fall (1) used conventional pistol bullets, which have an elliptical ogive with a length/diameter 
ratio close to unity. This rounded profile meets the symmetry axis with zero slope, so it is 
similar to the spherical profile modeled here (Figure 3). The pellets used by Kalista et al. (2) had 
somewhat blunted conical points, whereas the point of rifle bullets is often a spitzer ogive, which 
is a surface of revolution of a circular arc that is tangent to the shank. Figure 4 is an example of 
this kind of point, although in this case, the cylindrical shank of the bullet was not modeled. 
Finally, Fall showed that flat-ended cylinders left a hole that did not heal completely. A 
simulation of this case is shown in Figure 5. The response in the vicinity of the sharp corner of 


2 



this projectile is not captured correctly; attempts to use a finer mesh seemed to lead to numerical 
problems and the calculation could not be completed. 

As mentioned in the preceding section, the plate deformation in Figures 3-5 is confined to a 
fairly small region (not much beyond the projectile radius). The region of the very highest 
strains is a thin layer (less than 1 mm) immediately adjacent to the projectile. Figure 3 gives a 
good indication of the role of ductility in the penetration process. Even though the “bullet” has 
traveled a distance equal to nearly twice the initial thickness of the plate, the maximum principal 
strain (which occurs under the sphere) is still less than the tensile failure strain of Surlyn. In 
Figure 5 on the other hand, the thinnest region of the plate is adjacent to the corner of the 
cylinder. It is easy to imagine that continued thinning would lead to detachment of a plug of 
material. This failure mode would explain the lack of healing in this case (which is what Fall 
hypothesized was occurring with blunt cylindrical projectiles). 

Most published ballistics work involves metals and concrete. Perforation mechanisms have been 
reviewed (7) and include fragmentation, spalling, petalling, and plugging. A few studies on 
penetration of polymethylmethacrylate (8,9) and polycarbonate (10) indicate that the failure 
geometry depends on projectile shape, panel thickness, and strain rate. After being struck by a 
pistol bullet, some healed Surlyn plates contain fine radial lines that suggest that petalling may 
have occurred. This failure mode cannot, of course, be captured by a 2D simulation, and should 
be a subject for future work. In addition, a more sophisticated description of the post-yield 
mechanical behavior will probably be needed eventually: the constitutive model used in the 
present report includes neither strain-hardening nor damage accumulation. A final complication 
not treated here arises from the likelihood that the texture of semicrystalline polymers may 
evolve differently under multiaxial (as opposed to uniaxial) stretching. This means that material 
properties obtained in simple tension may not be relevant to the puncture experiment (11). 

Energy Absorption: Many approximate treatments of terminal ballistics have been proposed. 
Most closed-form expressions are either semiempirical energy-balance fits or they are based on 
the solution for the stress at the surface of an expanding spherical cavity (12). To employ the 
latter, the radial cavity stress is decomposed into tangential and axial components. Integrating 
the axial stress over the nose of a hemispherical projectile gives the assumed force on the nose: 

F = jta 2 {A+p t [(2/3)a(dV/dt) + 3V 2 /4]} (1) 

Where a is the hemisphere radius, V is the rigid-body projectile velocity, t is time, and p t is the 
target density (14). The elastic-plastic material constant A = (2a ys /3)[l + In (2E/3a y J. 
Neglecting friction and using the properties in Table 1 and an initial velocity of 260 m/sec, one 
calculates a deceleration of 2.9x1 0 2 m/sec/m. For a plate whose thickness is comparable to the 
bullet diameter (as in the experiments in references 1 and 2), the average deceleration would be 
expected to be lower, for the following reasons. In the early stages of penetration, the contact 
area would be less than the full hemisphere that was used in the closed-form solution, and in the 
later stages, the plastic region would not be constrained by elastic material (which is the implicit 
assumption of the spherical cavity solution). 
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Next, deceleration of a bullet was simulated using finite elements. A spherical projectile was 
given an initial velocity and mass that match those of the copper-jacketed 9 mm slugs used in 
experiments at LaRC (13). Figure 6 shows simulated projectile velocity as a function of distance 
for three different panel thicknesses calculated using the polymer properties in Table 1. For the 
thickest panel, the initial slope corresponds to about l.lxlO 3 m/sec velocity change per meter of 
travel. This is lower than the closed form result, as expected, and has the right order of 
magnitude. The sensitivity of this FE result to the material properties was explored using the 
variations in Table 3. These results are shown in Figure 7. 


Table 3. Constitutive parameters for sensitivity study in Fi gure 7. 



Modulus, 

MPa 

Yield Stress, 
MPa 

Case I 

350 

15 

Case II 

700 

15 

Case III 

350 

30 

Case IV 

700 

30 


The curves in Figures 6 and 7 are fairly featureless, except for a decrease in slope as the 
projectile nears the back surface of the plate and the polymer thins out. The sensitivity study 
shows that increases in either the polymer’s modulus or its yield stress would lead to more rapid 
deceleration of the projectile. These material properties will be revisited in the next section. 

The elastic-plastic model used here does not include a failure criterion, so it is not possible to 
predict the total expected velocity change with any accuracy. Upon examining Figure 3, 
however, it seems unlikely that a 6-mm panel would provide much additional resistance beyond, 
say, 12 mm of projectile travel. This would imply, from the FE results, a final velocity on the 
order of 250 m/sec. 

Experimentally, the velocity changes are much larger (Table 2). It therefore seems that factors 
omitted from the simulation — friction, viscoelastic dissipation, elastic dishing of the plate, and 
possibly boundary effects — play large roles in slowing the projectile. 


Table 2 Slowing of 9-mm Pistol Bullets by Surlyn (Ref. 13) (averages of 2 trials) 


Panel thickness, 
mm 

Initial velocity, 
m/sec 

Final velocity, 
m/sec 

6.477 

249 

203 

4.55 

263 

225 

2.26 

258 

239 


Rate and Temperature Effects: The time- and temperature- varying behavior of amorphous 
polymers near Tg has received a great deal of study. Consider first the time (rate) dependence. 
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The model proposed by Matsuoka (15) suggests that glassy modulus should scale with a power 
of the strain rate. The exponent should be about m/2, where m » 1/3 is the power law exponent 
of a creep curve. Although Surlyn is not glassy at room temperature, the creep compliance data 
of Fall (1) are consistent with m = 0.39. The nominal strain rate of the ballistic penetration is 
about 30xl0 3 sec' 1 , whereas the ASTM D638 tensile test (source of the data in Table 1) occurs at 
~1 sec' 1 . The effective modulus at ballistic rates could therefore be (30xl0 3 ) 1/6 =5.6 times higher. 
Low-rate puncture tests conducted on Surlyn at various rates (3) imply that the yield stress a ys is 
also likely to be at least twofold higher at the higher rate. For comparison, a similar factor of 2 
would be inferred from rate-dependent tensile yield data on polycarbonate (16). 

When plastic work is done on a material, a fraction of that work (often symbolized by (3) is 
converted to heat. For metals (3 » 0.9, but for polymers it tends to be lower (17) and (3 « 0.5 is an 
appropriate rule-of-thumb. The storage of “plastic” (actually anelastic) strain energy is the 
reason that polymer strains of at least 60% can be nearly completely recoverable (18). 
Representative values were chosen for polymer thermal conductivity (0.172 W/mK) and heat 
capacity (1.6 kJ/kgK). Incorporating heat generation with (3 = 0.5 in the FE calculations, it was 
found that temperatures climbed more than 50°C above ambient in a narrow zone near the 
projectile (Figure 8). This prediction matches reasonably well with the temperature rise seen 
experimentally (1). 

The properties of polymers are strongly dependent on temperature, so this temperature rise 
should, in turn, affect the mechanical properties. For this parametric study, the high-rate 
modulus and yield stress derived above were both allowed to decline linearly with temperature, 
reaching zero at the polymer melting point (~90°C). Figure 9 shows the calculated effect of 
thermomechanical coupling on energy absorption. 

Conclusions 


It was shown that polymer deformation under projectile impact is highly localized and should 
produce very localized heating. With a rounded projectile, strains are highest below the nose of 
the bullet and a highly ductile polymer may exhibit different failure modes than metals do. 
Directions for future work have been identified. 
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Contour Lines of materi strain t, Si-materi strain t. 

Deformation ( xl): mesh deform of TIME ANALYSIS, step 0.0025 



k 


x 


m 0 29951 
= 0.26623 

— 0.23295 
B 0.19967 
_ 0.16639 
S 0.13311 

0.099835 
1 0.066557 
M 0.033278 

— 0 



Figure 1. Low-rate impact (2.6 m/sec) on polymer plate. Sphere is moving downward. 


Contour Lines of materi strain t. Si-materi strain t. 

Deformation ( xl): mesh deform of TIME ANALYSIS, step 2.4e-5. 


k 




Figure 2. High-rate impact (260 m/sec) on polymer plate. 
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Contour Lines of materi strain t, Si-materi strain t. 

Deformation ( xl): mesh deform of TIME ANALYSIS, step 4e-5. 


k 


x 



m 4.7361 
M 4.2099 
B 3.6836 
5 3.1574 
E 2.631 2 
2.1049 
1.5787 
| 1.0525 
H 0.52623 


Figure 3. Deformation by impact of sphere on elastic-plastic plate. 


Contour Lines of materi strain t, Si-materi strain t. 

Deformation ( xl): mesh deform of TIME ANALYSIS, step 3.3e-5. 


kx 

Figure 4. 



_ 6.8736 
M 6.1098 
§j 5.3461 
4.5824 
g 3.8186 
: 3.0549 
2.2912 
g 1.5275 
M 0.76373 
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Contour Lines of materi strain t, Si-materi strain t. 

Deformation ( xl): mesh deform of TIME ANALYSIS, step 3.6e-5. 


k 


X 




Figure 5. As in Figures 3 and 4, but with cylindrical projectile. 



Figure 6. Calculated velocity of spherical projectile vs. distance for three plate thicknesses. 
Solid line: 6 mm plate; dotted line, 4 mm; dashed line, 2 mm. 
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Figure 7. Effect of polymer properties listed in Table 3 on deceleration of a spherical projectile 
impacting a 6-mm-thick plate. Solid: Case I; dotted, case II; short dashes, case III; long dashes, 
case 


Contour Lines of temp. 

Deformation ( xl): mesh deform of TIME ANALYSIS, step 4.5e-S. 



IV 


k 


x 


_ 56.216 
g 49.97 
3 43.724 
1 37.478 
_ 31.231 

— 24.985 
18.739 

— 12.493 
_ 6.2463 

^ -1 4489e-09 



Figure 8. Temperature contours (°K above ambient) for impact. Temperature- independent 
material properties. 
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Figure 9. Effect of thermomechanical coupling. Solid line: constant modulus and yield strength; 
dotted line: with softening due to plasticity-induced heating as described in the text. 
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Appendix 


A listing of a TOCHNOG script follows. Comments in parentheses are ignored by the program. 

( axispherelO ) 

(initiation section) 

(echo script to screen) 
echo -yes 

number_of_space_dimensions 2 

(declare variables to store in each node) 

( condif ^conduction and diffusion of heat) 

condi f_temperature 

materi_velocity 

mater i_s tress 

(plastic strain) 

materi_strain_plasti 

materi_strain_total 

( kappa=ef f ective accumulated scalar plastic strain) 
materi_plasti_kappa 

end_initia 

( bullet ) 

(point centered on circle with 4.51mm tolerance) 

geometry_point 1 0 4.5e-3 4.51e-3 

(mesh a circle and refine mesh twice) 

control_mesh_macro 1 -circle 0 2 

(set boundary nodes for detecting contact) 

control_mesh_macro_set_node_boundary 1 —yes 

(locate circle at [0,4.5mm] with 4.5 mm radius) 

control_mesh_macro_parameters 1 0 4.5e-3 4.5e-3 

(location to follow bullet movement during postprocessing) 

post_point 1 0 4.5e-3 

(bullet material properties: stiff, elastic) 
group_type 0 -materi 
group_materi_elasti_young 0 lei 3 
group_materi_density 0 19.58e3 
group_materi_elasti_poisson 0 0.3 
group_axisymmetric 0 -yes 

(constrain x-coordinates of bullet nodes (zero velocity)) 
bounda_unknown 13 -geometry_point 1 -velx 
bounda_time 13 0 

(apply initial velocity in y-direction) 

(select nodes to apply to) 
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control_data_put 15 -node_dof -geometry_point 1 

(assign initial values to all variables) 
control_data_put_double 15 0. (temp) 

0. -260 (Vx and Vy ) 


0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

( stress ) 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

( plas strain ) 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

( totl strain) 


0 . ( kappa ) 

(plate) 

group_type 1 -materi -condif 
group_condif_density 1 950 
group_condif_conductivity 1 0.172 
group_condif_capacity 1 1.6e3 

group_materi_elasti_young 1 1960e6 

group_materi_plasti_vonmises 1 30e6 
group_materi_density 1 950 
group_materi_elasti_poisson 1 0.4 
group_materi_plasti_heatgeneration 1 0.5 
group_materi_memory 1 —updated 

(Young's modulus and yield stress dependence on temperature) 
dependency_item 4 -group_materi_elasti_young 1 -temp 3 
(piecewise linear: from temp=0 (ambient) to 70; 70 to 200) 
dependency_diagram 4 0 70 200 1960e6 0 0 

dependency_item 6 -group_materi_plasti_vonmises 1 -temp 3 
dependency_diagram 6 0 70 200 30e6 0 0 

group_axisymmetric 1 -yes 

control_mesh_macro 20 -rectangle 1 10 4 
(plate is 15mmx6mm) 

control_mesh_macro_parameters 20 7.5e-3 -3.1e-3 15e-3 6e-3 

control_mesh_macro_element 20 -quad4 
control_mesh_macro_set_node_boundary 20 -yes 
(clamped edge) 

geometry_line 3 15e-3 -10e-3 15e-3 10e-3 0.1e-3 
(center line) 

geometry_line 4 0 -7e-3 0 -0.1e-3 0.1e-3 
(boundary conditions) 

bounda_unknown 5 -geometry_line 3 -vely -velx 
bounda_time 5 0.0 0.0 100.00 0.0 

bounda_unknown 7 -geometry_line 4 -velx 

bounda_time 7 0.0 0.0 100.00 0.0 

contact_penalty_velocity le7 

(turn off "stabilization" (useful for quasistatic case only) ) 
options_stabilization -no 
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(generate completely new mesh for plate; 0.5 mm elements) 
control_mesh_new_mesh 30 .0005 
control_mesh_new_mesh_element 30 -quad4 
control_mesh_new_mesh_region 30 1 

(take 50 steps of 0.1 microsecond) 
control_timestep 40 le-7 50e-7 
(explicit calculation) 
control_timestep_iterations 40 1 

(print time and projectile speed to screen for each step) 
print_filter 1 -post_point_dof 1 -vely 
control_print 40 -time_current -post_point_dof 
(save data to file) 

control_print_data_versus_data 40 -time_current 0 0- 
post_point_dof 1 2 

(save output files for visualization) 
control_print_gid 45 -separate 

(loop back to statement 30 and repeat 8 more times) 
control_repeat 50 8 30 

end data 
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